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ABSTRACT 

The detection of multiple mode harmonic kink oscillations in coronal loops enables to obtain information 
on coronal density stratification and magnetic field expansion using seismology inversion techniques. The 
inference is based on the measurement of the period ratio between the fundamental mode and the first overtone 
and theoretical results for the period ratio under the hypotheses of coronal density stratification and magnetic 
field expansion of the wave guide. We present a Bayesian analysis of multiple mode harmonic oscillations for 
the inversion of the density scale height and magnetic flux tube expansion, under each of the hypotheses. Then, 
the two models are compared using a Bayesian model comparison scheme to assess how plausible each one is, 
given our current state of knowledge. 

Subject headings: magnetohydrodynamics (MHD) — methods: statistical — Sun: corona — Sun: oscillations 
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1. INTRODUCTION 

Observations of magnetohydrodynamic (MHD) oscil- 
lations in solar coronal loops indicate the simulta- 
neous presenc e of multiple mode harmonics in the 
same structure ( Verwichte et al.' '2004; De Moortel & Brady] 
2^ i Van D oorsselaere et al. 2007; O'Sheaetal. 2003 
van Doorsselaere et al. 2009). Their potential use for the di- 
agnostic of the plasma co nditions in the coron a has attracted 
considerable interest (see lAndries et al.ll2009[ fo r a rev iew). 
The idea was first put forward by lAndries et al.l (120051) and 
lOoossens et al. ( 2006) who found that the ratio of the fun- 
damental mode period to twice that of its first overtone in 
the longitudinal direction depends on the density structur- 
ing along magnetic field lines. It is equal to unity in longi- 
tudinally uniform and unbounded tubes, but is smaller than 
one when density stratification is pres ent. Using o b serva- 
tional measurements of period ratios, Andries et al.l (12005b 
andlVan Doorsselaere et al. (2007) obtained estimates for the 
coronal density scale height. 

An alternat ive hypothesis was formulated by 

IVerth & Erdelvil (120081) whereby the expansion of the 
magnetic loop produces a deviation from unity in the period 
ratio, increasing its value. The effect is relevant for obtaining 
accurate estimates o f the density sca le height using multiple 
mode period rat i os (iVerth & Erdelyii2008; , Verth et al.,2008i: 
iRuderman et al.l l2008h . Observational evidence for mag- 
netic t ube expansion has been reported by e.g., iKlimchuM 
(120001) in soft X-ra y loops observed with Yohkoh and by 
IWatko & Khmc huk (2000) in non-flare and postflare loops 
observed with TRACE. Observati ons of period r atios l arger 
than one have been report ed by lO'Shea et al.l (120071) and 
iDe Moortel & Bradyl (l2007l) . 

Since then, period ratio studies have analyzed differ- 
ent models for the density structuring in coronal loops 
(jMcEwan et al. 2006, 2008), the influence of the elliptic shape 
of loops (Morton & Erdelvi 2009), the twist of the mag- 
netic field (Karami & Baharii2012) . the effect of the environ- 
ment (lOrza et al.ll20T2ir and the temporal expansion of loops 
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dBallai & Orzall20T2h . Period ratios have also been a nalyzed 
for slow MHD modes (Macnamara & Roberts 2010) and for 
sausage modes (.Macna mara & Roberts 20 11) and their use 
has been suggested in the context of prominence seismology 
dDiaz et alj|2010l; lArregui et al.ll20T2l) . 

The appUcation of Bayesian analysis techniques to coro- 
nal seismology is in its infancy. A first attempt of parameter 
inference using damped loop osc illations was presented by 
lArregui & Asensio RamosI (1201 ll) . In this paper, we present 
the first application of Bayesian model comparison techniques 
to coronal seismology. We first perform Bayesian parameter 
inference for the coronal density scale height and the mag- 
netic tube expansion, under the two hypotheses of density 
stratification and magnetic field divergence, using multiple 
mode oscillations. Then, we assess which one of the two hy- 
potheses better explains the observations, for given values of 
the period ratio. 

2. THEORETICAL MODELS 

The deviation from unity for the ratio between the funda- 
mental and the first overtone transverse kink oscillation peri- 
ods has been attributed to two main physical effects. Each one 
constitutes a hypo thesis to expla i n the d ata. 

In the model by lAndries et al.l (120051) coronal density strat- 
ification produces a decrease of the period ratio. This model 
projects a vertically stratified isothermal atmosphere onto a 
semicircular loop. An analytical expression for the depen- 
de nce of the pe riod ratio on density scale height was obtained 
bv lSafari et afl (12007) . This expression can be rewritten as 
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with Pi and P2 the periods of the fundamental and first over- 
tone modes, r\ - L/nH the ratio of the loop height at the apex 
to the density scale height H, and L the loop length. We have 
checked that Equation ([T]) provid es us with a good ap proxi- 
mation to the num erical results by lAnd ries et al.l (120051) . 

In the model by IVerth & Erd elvi (2008]) magnetic tube ex- 
pansion produces an increase of the period ratio given by 



= 1 + 



3(F2 - 1) 



(2) 



2P2 ' ' 27r2 ' 
with the expansion defined as F - ra/r/, where is the radius 
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Figure 1. (a) Posterior distributions for rj, under tlie density stratification model, for two values of the period ratio, (b) Posterior distributions for the magnetic 
tube expansion, T, under the magnetic expansion model, for two values of the period ratio. The measured period ratio and the inferred median of the distribution, 



: 0.91 ±0.04, = 1.26+ 
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1.87+° °?. The improved errors on r in the measurements of lVerwichte et al.l <2004l) have been taken from lVan Doorsselaere et all i2007l) . 



with uncertainties given at the 68% credible interval are: (a) r 

r= 1.38 ± 0.04, r 

Similar errors are assigned to the mean values in the measurements bv lDe Moortel & Bradvl <2007l) . 

at the apex and rj is the radius at the footpoint. Equation ^ 
was obtained under the assumption that |r - 1 1 1 . More- 
over, it is the result for a particular expanding loop model. 
Another equilibriu m state will produce a qu antitatively differ- 
ent result, although iRuderman et alJ (120081) anticipate a qual- 
itatively similar result. 

Both longitudinal stratification and magnetic tube expan- 
sion have forward models that relate one observed quantity, 
the period ratio P1/2P2, to one physical quantity to be in- 
ferred, 77 or r. 



0.79 ± 0.03, = 3.39+° '^. (b) r = 1.07 ± 0.04, r = 1.20^ 



3. BAYESIAN PARAMETER INFERENCE 

To perform the inference using existing estim ates for the 
period ratio we employ Bayes' theorem ( Baves & Pricell763l) 



p(e\D, M) = 



p{D\e,M)p{e\M) 
Jd0p(D\0, M)p{0\My 



(3) 



which gives the solution to the inverse problem in terms of the 
posterior probability distribution, p{6\D,M), that describes 
how probability is distributed among the possible values of 
the unknown parameter, 0, given the data D and the assumed 
model M. The function p(D\0, M) is the likelihood of obtain- 
ing a data realization actually observed as a function of the 
parameter vector and provides a measure of how well the data 
are predicted by the model. The prior probability p{0\M) en- 
codes any prior information we might have on the model pa- 
rameters, without taking into account the observed data. The 
denominator is the evidence, an integral of the likelihood over 
the prior distribution. This quantity plays no role in parameter 
inference, but will become central in the model comparison 
described in Section|4] 

3.1. Longitudinally Stratified Loops 

For longitudinally stratified loops, model Mi , the forward 
problem is given in Eq. ([T]l. To evaluate the likelihood, we 
assume the model is true. Then, the period ratio measurement 
(r) will differ from the prediction (ri) because of measure- 
ment uncertainties (e), so that r = ri + e. The probability of 
obtaining the measured value is equal to the probability of the 
error Assuming Gaussian errors, the likelihood for model Mi 
is then expressed in the following manner 



p(r\7]. Ml) 



V2^c 



■ exp 



(r-rif 



2cr2 



(4) 



with cr^ the variance associated to the observed period ratio. 
In the following we assign observed period ratio errors to the 
standard deviation cr. 

The prior indicates our level of knowledge (ignorance) be- 
fore considering the observed data. We have adopted a uni- 
form prior distribution for the unknown, 77, over a given range, 
so that we can write 



pinlMi) = 



1 



for 1]" 



<n < 



(5) 



and zero otherwise. We only consider loops with density de- 
creasing with height. We know from observations that very 
large coronal loops with the apex height above several atmo- 
spheric scale heights exist, so we consider the range 77 e [0, 8] 
in the prior above. 

Parameter inference is next performed using Bayes' theo- 
rem (l3]l. We use period ratio measu rements reported in ob- 
servations by ' Verwic hte et all (120041) . We must note that un- 
ambiguous identification of measured periods with particular 
overtones is a difficult task. Early observations were limited 
to the measurement of different periodicities at a single point, 
without spatial information ( Verwichte et al. 2004). Even in 
the case in which information at differen t cuts along the loops 
is available dPe Moortel & Bradvl2007l) . it is unclear how one 
should assign those periods to theoretical overtones. As our 
paper aims at providing a method to perform parameter infer- 
ence and model comparison, we have considered peri od ratio 
values discussed in the review bvl Andrieset"an(l2009h . notmg 
that the reliability of the results is closely related to that of the 
adopted assumptions and theoretical interpretations. 

Figure [T^ shows posterior probability distributions for 77 
computed usingEq. (|3]l with likelihood and prior given by 
Eqs. dUi and o}, for two period ratio measurements by 
IVerwichte et al.rii2004l) . Well constrained distributions are ob- 
tained. For the measured period ratios r ~ 0.79 and r ~ 0.91, 
the inversion leads to density scale heights of // = 21 Mm and 
// = 56 Mm, respectively, for a loop with a height at the apex 
ofL/jT = 70 Mm. 
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3.2. Expanding Magnetic Loops 

For expanding loops, model M2, the forward problem is 
given in Eq. (|2]l. To evaluate the likelihood, we follow the 
same procedure as above, which leads to 



p(r\r,M2)^ 



1 



■ exp 



2cr2 



(6) 



We also adopt a uniform prior distribution for the unknown 
r, so that we can write 



Pimi) 



1 



for r™" < r < r" 



(7) 



and zero otherwise. We consider F e [1, 2.5]. 

Figure \^ shows posterior probability distributions for F 
computed using Eq. Q with likelihood and prior given by 
Eqs. ^ and (|7]l, for two period ratio measure ments. Again, 
well c onstrained distributions are obtained. lAndries et al.l 
(120091) discuss period ratio measurements in Table 1 by 
iDe Moortel & Bradvl (12007b . Assuming that either the most 
power is in the fundamental mode or in the first overtone, 
mean values for the period ratio of r ~ 1.07 and r ~ 1.38 
are obtained. For those values, the inversion leads to tube 
expansion factor s that are co mpatible with the estimates by 
iKlimchukI ( 12000 ) and Watko & KlimchukI ( 120001) . Note how- 
ever that, according to Fig[TJ), a period ratio of r ~ 1.38 re- 
quires an expansion of the tube by a factor of F ~ 1.85, while 
observations by Wa tko & KlimchukI (12000) seem to indicate 
that in only very few cases does this parameter approach or 
exceed a value of 2. 

For both theoretical models, the Bayesian framework 
makes use of all the available information in a consistent man- 
ner and enables us to consistently propagate errors from ob- 
servations to inferred parameters. 

4. BAYESIAN MODEL COMPARISON 

Parameter inference was performed under the hypothesis 
that observed period ratios are the effect of either density strat- 
ification or magnetic tube expansion. The probability distri- 
butions in Figure[T]are conditional on the specific models used 
to explain the observations. The Bayesian framework enables 
us to present different models to the same data and assess 
in a quantitative manner which one is favored by them. We 
have compared three models: Mo for a uniform density and 
magnetic field strength tube. Mi for a longitudinally stratified 
loop, and M2 for a expanding magnetic loop. 

As model Mq predicts a period ratio of one, regardless of rj 
or F, we can write 



p(r\Mo) 



1 



^/2nc 



■ exp 



(r-iy 

2cr2 



(8) 



for the likelihood in this case. Note that p(r\MQ) - 
p(r\T],Mo)^ p(r\T,Mo). 

To determine the plausibility of models Mi and M2 between 
them and with respect to model Mo, we evaluate the posterior 
probabilities to ascertain the relative merits of two models, M, 
and Mj. This is done by applying Bayes' theorem (Eq.O to 
the two mod els and considering posterior ratios of the form 
(iJeffrevsll 19611) 



p(Mi\r) ^ p(r\Md p(Mi) 
p(Mj\r) p(r\Mj)p(Mjy 



(9) 



The first ratio on the right-hand side is the Bayes factor It 
expresses how well the observed data are predicted by model 
Mi, compared to model My. The second ratio, the prior odds 
ratio, measures how much our initial beliefs favored M, over 
Mj, before considering the data. As we have no particular a 
priori preference for one model over the other, before consi d- 
ering the data, we take piM,) - p{Mj) - 1/2. Our assessment 
of the plausibility of models is then based on the computation 
of the Bayes factor of M, against Mj given by 



BFi. 



p(r\Md 
piAMj) 



(10) 



In model selection, we are interested in the most proba- 
ble model, independently of the parameters, i.e., we should 
marginalise out all parameters. This is achieved by perform- 
ing an integral of the likelihood over the full parameter space. 
The marginal likelihood for a given model M, is then 



/7(r|M,)= pir,e\Mi)d9^ p(r\0, Mi)p(0\Mdde, 

Jgmin Jgrnii, 

(11) 

where e [0™'", g'™*] represents the particular parameter of 
the model and we have used the product rule to expand the 
probability of r and 0, given model M, . 

Figure |2^ displays the marginal likelihoods for the three 
considered models. For a given observed period ratio, the 
plausibility of one model over the other is given by the ra- 
tio of these two quantities at the measured period ratio r. The 
uniform model has a marginal likelihood that is maximum at 
one. The model with density stratification is clearly favored 
for observations of period ratios below unity. Figure |2j) quan- 
tifies the relative performance of models M\ and Mo by com- 
puting the Bayes factor The lower r, the more plausible Mi 
is aga inst Mq. Jeffreys' scale (Jeffreys 1961; Kass & Rafteryl 
assigns different levels of evidence to the values of the 
Bayes factor in natural logarithm units. They are shown us- 
ing different darkening options in Figures |2j>-d. According to 
Figure |2j3, and given the assumed uncertainty of cr = 0.08, 
a period ratio measurement should be considered as positive 
evidence for model M\ against model Mq only if it is below 
0.87. As strong evidence, only if it is below 0.78. A period 
ratio below 0.71 indicates very strong evidence for model Mi . 

A similar comparison was made for models M2 and Mo 
(Figure|2};). The model for magnetic tube expansion is clearly 
favored for observations of period ratios above unity. Model 
Mo cannot be ruled out by the simple fact that r > 1 . Positive 
evidence for model M2 against model Mo exists for period 
ratio measurements above 1.16. If r > 1.23, there is strong 
evidence for model M2. A period ratio above 1.28 indicates 
very strong evidence for model M2. 

Finally, we have compared model Mi against model M2, in 
view of observed data. Figure |2}l shows that the evidence is 
the same for both models if the period ratio is unity. Around 
this value, there is a region, 0.96 < r < 1.06, in which no 
firm conclusion can be established. Positive, strong, and very 
strong evidence for Mi occur below 0.96, 0.87, and 0.80, re- 
spectively. Positive, strong, and very strong evidence for M2 
occur above 1.06, 1.15, and 1.21, respectively. 

In our model comparison, cr - 0.08 has been selected so 
as to clearly show the different regimes for the evidence. An 
increase (decrease) of cr produces a decrease (increase) of any 
evidence. The inference in Figure[T^ for r - 0.91 with uncer- 
tainty of cr = 0.04 falls into the region of positive evidence for 




Figure 2. (a) Mai'ginal likelihoods computed using Eq. jl U for models Mq, M\ and M2 as a function of the data, (b)-(d) Bayes factors computed using Eq. UOt 
for the model comparisons between M\ over Mo, Mi over Mo, and M\ over Mt. In (b)-(d), white regions indicate evidence not worth more than a bare mention 
(2 logpBF e [0, 2]); PE: positive evidence (2 iog^BF e [2, 6]); SE: strong evidence (2 log^BF e [6, 10]);VSE: very strong evidence (2 loggBF > 10). Uniform 
priors in the ranges t] e [0, 8] and E e [1, 2.5] have been taken. In all figures cr = 0.08. 



hypothesis Mi . The inference in FigurefTJi for r - 1 .07 should 
have an uncertainty of cr = 0.03 (close to the reported error) 
to be considered done under positive evidence for hypothesis 
M2. 

5. CONCLUSIONS 

We have presented a simple and straightforward technique 
to perform parameter inference and model comparison us- 
ing period ratios of kink oscillations in the Bayesian frame- 
work. Parameter inference enable us to obtain estimates for 
the coronal density scale height and the magnetic loop expan- 
sion, using all the available information and with correctly 
propagated uncertainties. Density scale height estimates are 
compatible with previous studies. Magnetic tube expansion 
fact ors are compatible w i th the estimates by Klimchuk (2000) 
and lWatko & KUmchukI ( 12003) . Even if Eqs. ([B and ^ im- 
ply well-posed inversion problems from a mathematical point 
of view, the unknown parameters cannot be obtained uniquely, 
due to the uncertainties associated to observed data. The 
Bayesian framework enables us to consistently deal with this 
problem. 

Bayesian model comparison enable us to assess the plausi- 
bility of the inferences, which are conditional on theoretical 
models. Deviations of the period ratio below (above) unity do 
not necessarily imply density stratification (tube expansion) 
to be preferred in front of the null hypothesis. The degree of 
evidence for the two theories so far invoked can be quanti- 
tatively assessed with the use of Bayes factors, to precisely 
decide which one of the two hypotheses is more plausible, for 



a given measured period ratio and the associated uncertainty. 

Our analysis provides a simple and easy to use method to 
perform inference and model comparison in the presence of 
incomplete and uncertain information. Measured period ra- 
tios and their uncertainties determine the strength of the evi- 
dence in favour of a particular hypothesis and, therefore, the 
support of the performed inferences. 

The Bayesian formalism for inference and model com- 
parison is the only fully correct way we have to get in- 
formation about physical parameters and the plausibility 
of hypotheses from observations (see e.g. iTrottal I2008t 
I von Tou ssaind 1201 ll) under incomplete a nd uncertain infor- 
mation. Arregui & A sensio RamosI (1201 ll) performed the first 
Bayesian parameter inference in coronal seismology. This 
paper presents the first application of Bayesian techniques 
to model comparison. When analytical forward problems 
are available, both parameter inference and model compari- 
son reduce to a simple and straightforward evaluation of the 
marginal posteriors and the Bayes factors, that can be ob- 
tained by solving simple integrals. The methods outlined in 
this paper can be directly applied to most of the seismology 
inversion problems in which other physical effects and param- 
eters are involved. For instance, additional effects that influ- 
ence period ratios could be compared to the already consid- 
ered hypotheses. 
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